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INTRODUCTION 


It  has  long  been  known  that  a  shock  wave  is  propagated  in  a  gas 
with  a  velocity  which  changes  with  pressure  and  which  nay  be  many  times 
greater  than  the  velocity  of  sound.  Accordingly,  the  velocity  is 
measured  by  the  pressure  and  the  energy  content  in  the  wave,  and  is, 
therefore,  a  rather  significant  factor  in  the  determination  of  the 
demolition  capacity  of  the  wave.  Nevertheless,  knowledge  regarding 
shock  velocity  is  quite  incomplete  and  unsatisfactory  despite  the  many 
investigations  carried  out  in  the  past. 

The  measurements  presented  here  "-ere  made  at  the  Fysikaliska 
Forskningsa^elningen  vid  AB  Bofors  [physical  Research  Department  at 
Bofors,  LtdJJ  ,  and  aim  to  cope  direcuiy  with  this  deficiency  regarding 
spherical  charges. 

I-Iuowledge  of  the  results  also  yields  the  opportunity  to  check 
theory  by  examining  the  connexion  between  the  emanating  wave  and  the 
characteristics  of  the  explosive  cnarge.  This  was  another  reason  for 
the  present  rather  comprehensive  investigation. 

'<7e  7/ere  particularly  anxious  to  check  two  theories:  an  old  or.e 
by  Rudenberg  (Reference  l),  and  a  quite  recent  one  by  Sedov  (Reference 
2)  because  of  application  to  the  calculation  of  the  amount  of  energy 
that  is  liberated  by  an  atomic  bomb  explosion.  Belov/  are  shown  the 
Bofors  measurements  which  are  not  at  all  in  agreement  with  the  pre¬ 
dictions  of  either  of  the  above  theories.  Thus  it  is  certain  that  the 
above  theories  are  useless  for  the  calculation  of  the  amount  of  energy 
liberated  by  an  atomic  bomb  explosion.  Applications  to  the  atomic  bomb 
appear,  therefore,  to  be  hardly  promising. 

Let  us  survey  the  previous  experimental  work  published  on  this 
topic.  One  then  meets  only  with  qualitative  results.  Many  investi¬ 
gators  have  found  that  shock  front  velocity  is  very  high  in  the 
vicinity  of  the  charge  and  that  it  decreases  with  distance,  gradually 
approaching  acoustic  velocity,  and  from  their  findings  have  attempted 
to  formulate  laws,  but  the  measurements  turned  out  to  be  irreproducible . 

The  essential  cause  of  the  poorness  of  these  results,  despite  the 
range  of  weights  of  charges  exploded,  which  in  some  cases  extended  up  to 
1500  kg,  is  indubitably  due  to  the  neglect  of  the  effect  of  the  shape  of 
the  charge.  Measurements  at  Bofors  have  shown  that  cylindrical  charges 
produce  a  distinct  directional  effect  not  only  with  respect  to  the 
impulse,  as  the  author  has  demonstrated  in  Reference  3,  out  also  with 
regard  to  the  velocity.  Thus  it  is  apparent  that  when  velocity 
measurements  are  obtained  without  consideration  of  the  shape  of  tne 
charge,  the  results  will  be  rather  puzzling  and  irreproducible.  This 
remark  applies  to  all  the  experiments  carried  out  by  the  authors  named 
below.  It  should  further  be  added  that  the  charges  were  often  inserted 
in  magazine  containers  or  metallic  covers  which  could  have  directed  the 
explosion  wave  in  an  imperspicuous  manner. 


The  first  Telocity  measurements  of  explosion  waves  were  made  by 
Mach,  who,  in  1877  to  1889,  experimented  with  very  small  charge  weights. 
Among  other  experiments  we  may  list  those  by  Wolff  (Reference  4),  von 
Angerer  and  L&denburg  (Reference  5),  Burlot  (Reference  3),  and  Partlo 
and  Service  (Reference  7). 


The  two  last  named  experimenters  have  made  the  notable  observation 
that  the  velocity  within  a  certain  range  of  charge  radi:'  is  subsonic. 
It  does  not  appear  likely  that  this  result  is  correct.  Apparently,  the 
error  lies  in  the  neglect  of  the  directional  effect  and  of  the  scaling 
law  derived  from  the  principle  of  dynamic  similarity  according  to  which 
all  velocities  are  equal  at  distances  measured  in  charge  radii,  for  a 
specified  type  of  explosive. 

The  form  of  the  charge  introduces  considerable  complications  as  to 
the  directional  effect,  particularly  near  the  ground.  It  is  reasonable 
to  begin  by  eliminating  these  difficulties  by  choosing  as  a  charge  form 
the  sphere  which  alone  is  fully  symmetric  if  set  off  at  the  charge 
centre.  This  is  the  reason  for  confining  the  present  investigation  to 
spherical  charges. 


When  such  a  charge  is  detonated  it  3ends  out  a  shock  wave  within 
which  the  pressure  behaves,  as  Figure  1  shows,  with  a  pressure  maximum 
at  the  discontinuous  front.  The  pressure  generally  decreases  linearly 
with  time  behind  the  front  and  attains  at  point  2  the  value  of  one 
atmosphere  to  drop  below  this  value  thereafter.  The  wave  front  is 
propagated  with  a  velocity  w^  which  decreases  with  the  distance  from  the 
charge,  starting  from  the  maximum  value  wQ  at  the  surface. 


Rq  =  charge  radius 

R  =  distance  from  the  centre  of  the  charge 

T^  =  time  for  shock  front  to  reach  R,  starting  from  Rq 

T^  =  time  for  pressure  to  become  atmospheric  after  passage 
of  the  shock  front 

Q  =  charge  7/eight 

=  shock  velocity  at  time  T^. 


Then  the  corresponding  quantities  scaled^according  to  the  law  of  simili¬ 
tude  are: 


r  =R  Q 


ti  =  TXQ 


-V3 

i  y 

r  =  R  Q-1/3 

n 

-1/3 

>  y 

t2  =  T2q- ^ 

The  times  t,  and  t„  are  related  to  r  and  r  by 
12  o  J 

t«  =f,(r  -  r  ) 


ti=fi(r-ro),  u2  "*■  2  o' 


where  f^  and  f^  depend  on  the  type  of  explosive  but  do  not  depend 
explicitly  on  Q. 
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FIG.  I 
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Pressure-Time  Curve  for  Spherical  TNT  Charges 

FIG.  2 


Typical  Cathode-Ray  Oscillogram  Showing  Measured  Time  Intervals 

FIG.  3 
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The  difference  t+  =  t  -  t^  is  the  duration  of  the  positive  pressure 

at  a  fixed  point  r  while  a  wave  is  passing.  The  positive  wave-length  1 
is  the  distance  between  the  position  of  the  front  at  time  tg  and  r.  One 
then  gets 

T  W1  +  w2 

\  =w2(t2  -  t1),  [sic;  1+  = - 2 - (t2  -  tx)  is  a  better 


approximation^ 


where  w2  may  approximately  be  replaced  by  the  sonic  velocity  c  at  large 
distances. 


EXPERIMENTAL  PROCEDURE 

At  the  moment  the  explosion  products  break  through  the  surface  of 
the  charge,  almost  immediately  a  shock  wave  forms  and  an  intensive 
radiation  is  emitted  and  is  recorded  on  an  oscillograph  by  means  of  a 
highly  evacuated  photocell  (Philips,  type  3250)  and  a  trigger  coupling. 
At  the  distance  R  is  placed  a  crystal-receiver  which  registers  the  time 
of  arrival  of  the  wave  at  R. 

The  arrangements  have  to  be  varied  according  to  whether  the  dis¬ 
tance  is  small  or  large  because  in  the  former  case  one  measures  time  of 
the  order  of  microseconds  and  excess  pressures  that  are  exceedingly  high, 
while  in  the  latter  case  one  measures  longer  time  intervals  of  the.  order 
of  milliseconds  and  smaller  excess  pressures. 

At  small  distances  one  uses  an  S.E.  quartz-gauge.  It  consists  of 
a  steel  rod,  at  the  farthest  end  of  which  one  places  a  quartz  crystal 
shielded  by  a  threaded  hood.  The  oscillograph  is  a  polar-recording 
cathode-ray  "oscillograph  built  by  the  Fysikaliska  Forskingsavdelningen, 
A3  Bofors  which  contains  a  cathode-ray  tube  manufactured  by  the  A3 
Standard  Radio  Fabrik,  Stockholm.  Spots  of  light  describe  a  spiral 
consisting  of  about  three  effective  turns  with  a  speed  that  is  adjust¬ 
able  in  steps  from  20.000  to  200  rps.  For  a  shorter  time  interval  a 
higher  speed  of  rotation  of  the  oscillograph  trace  is  used.  Thus  cer¬ 
tain  difficulties  arise  on  how  to  start  the  oscillograph  at  the  right 
moment  so  that  one  loses  the  incoming  signals  from  neither  the  photo¬ 
cell  nor  the  quartz-gauge.  These  signals  are,  as  a  rule,  three  in 
number,  taken  at  various  distances  from  the  charge.  Thus  three  travel 
times  are  measured  for  each  explosion.  The  problem  is  solved  in  the 
following  manner :  The  charge  is  initiated  by  a  detonator  cap  and  a 
PSTN-fuze*  through  which  a  thin  copper  thread  is  drawn  out  to  a  suitable 
distance  from  the  explosive  cap.  TIhen  this  wire  breaks,  the  oscillo¬ 
graph  commences  to  sweep  its  spiral  path,  which  is  recorded  on  the 
screen.  Using  various  lengths  for  the  fuze  between  the  break-wire  and 
the  detonator  cap  in  the  centre  of  the  charge  one  may  easily  reproduce 
various  delay  times.  Such  an  oscillogram  is  shown  in  Figure  2.  The 
oscillograph  tubes  do  riot  indicate  full  circular  symmetry.  Therefore, 
the  reading  is  performed  with  the  aid  of  a  calibrated  protractor.  With 
a  speed  of  10.000  rps  the  reading  accuracy  is  about  0.5  jjus,  i.e.,  one 
may  measure  50  /u-s  with  a  tolerance  of  1  per  cent. 

•  [P5TH  =  05H3012S^ 
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Of  course,  it  is  important  for  such  short  tine  intervals  that  the 
response  tins  of  the  entire  circuit  be  minimal.  This  requirement  is 
checked  by  firing  several  shots  with  a  quartz  crystal  cemented  directly 
to  the  surface  of  the  charge  at'  each  firing.  The  signal  from  the 
crystal  and  the  photocell  trigger  circuit  are  incorporated  simultane¬ 
ously  in  the  cathode-ray  oscillograph  which  has  a  recording  velocity  of 
4000  m  sec.“^.  It  turned  out  that  the  time  difference  between  the  two 
signals  maintained  the  magnitude  of  1  to  2  times  10”^  sec.,  i.e.,  a 
value  safely  below  the  accuracy  of  reading. 

At  large  distances  a  six-spiral  Siemens  oscillograph  was  employed 
with  a  maximal  paper  speed  of  10  m  sec“^.  The  crystal  receiver  was  a 
Brush  microphone  with  Seignette  crystal.  The  microphone  is  connected 
to  a  non-linear  amplifier,  that  is,  one  with  very  high  sensitivity  for 
small  incut  voltage,  but  with  currer'4-  saturation  already  at  moderate 
voltage  [[i.e.,  a  high  gain  amplifier]  .  Accordingly,  zero  time  in  the 
pressure-time  curve  may  be  measured  with  an  accuracy  greater  than  the 
one  obtainable  with  the  use  of  a  linear  amplifier.  For  this  compare 
Figure  3. which  shows  the  distorted  and  normal  diagrams  according  to 
Figure  1.  The  current  curve  of  the  photocell  is  also  shown  in  the 
figure.  In  order  to  facilitate  distance  measurements  and  to  reduce 
effects  of  possible  horizontal  wind  disturbances  the  charge  is  suspended 
vertically  above  the  microphone  placed  about  1  m  above  the  ground.  Time 
measurements  are  made  with  a  tuning  fork  calibrated  with  cn°  from  the 
Telegraph  Company  at  a  telephonically  transmitted  frequency  of  500  Hz 
Q..e.,  cpsj,  the  accuracv  of  which  is  indicated  to  be  greater  than  1  in 
105.  “ 

Explosions  in  compressed  and  evacuated  air  are  carried  out  in  a 
specially  constructed  spherical  chamber  of  5.5  m^  displacement  — ith  a 
diameter  of  2  m.  This  permits  the  explosion  of  0.25  kg  of  TUT  H-  C^HgOgN^ 
at  an  initial  pressure  or  4  ata  Atmospheres  absolut^Q  .  The  corhr-iner  °  " 
can  also  be  evacuated  to  a  vacuum  of  a  few  mm  Eg.  The  measuring  arrange¬ 
ments  are  the  same  as  given  above  for  small  distances  in  free  air.  The 
photocells  are  located  outside  the  container  and  react  to  the  radiation 
from  within  through  a  small  glass  window.  Quartz-gauges  are  placed  in¬ 
side  the  chamber  and  are  connected  by  pressure-sealed  leads  to  each  of 
their  trigger  circuits  on  the  outside. 

EXPERIMENTAL  RESULTS  AND  CORRECTIONS 

Charges  are  weighed  in  grams,  small  distances  are  measured  in 
millimeters,  and  large  distances  in  centimeters.  Temperature,  humidity 
and  barometric  pressure  are  determined  before  every  shot.  The  charges 
are  exploded  in  calm  air. 

The  measurements  are  corrected  as  f ollows : 

1.  Quartz-gauge  measurements  of  travel  times  are  decreased 
by  1  ^s  corresponding  to  the  time  the  pressure  requires  to  be 
transmitted  through  the  protective  shield  of  the  gauge. 
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2.  Travel  tines  are  reduced  to  a  temperature  of  C  deg  C  and 
a  humidity  of  0  per  cent  by  multiplication  by  the  factor 

k  =  1  +  1.832  x  I0'  +  0.16  eh/3 

wher  e 

-S'  =  temperature  in  deg  C, 
e  =  humidity  in  per  cent, 
h  =  measured  vapour  pressure, 

B  =  barometric  pressure. 


The  factor  k  indicates  the  ratio  of  sonic  velocity  in  the 
given  still  atmosphere  to  that  in  an  atmosphere  of  0  deg  C  and  0  per  cent 
humidity. 


3.  The  times  measured  by  the  polar-recording  oscillograph  are 
reduced  by  0.39  per  cent  because  it  permits  agreement  with  the  more 
accurately  calibrated  time  indication  of  the  spiral  oscillograph. 

The  parameters  varied  within  the  following  limits: 


Temperature 

Humidity 


2.5  to  21.2  deg  C, 
61  to  96  per  cent. 


3arometric  Pressure  738  to  759  mm  Eg, 


Charge  freight 


Type  of  Explosive 


0.245  to  0.255  kg, 

0.965  to  1.005  kg. 

PSTN;  TUT  with  10,  20,  00  per  cent  Ai. 


The  corrected  mean  values  of  the  times  are  found  in  Tables 
I-toHEas  follows  s 


Table  I  gives  the  times  t^  and  t.,  at  atmospheric  pressure. 

Table  H  gives  the  times  t^  at  various  densities  of  the 
ambient  air. 

Table  TTT  gives  the  times  t^  at  atmospheric  pressure  for 
explosive  charge a  other  than  100  per  cent’ TUT, 
namely:  PET” }  TNT  with  1C,  20,  30  per  cent  Al. 

With  the  few  measurements  that  we  have  for  any  given  dis¬ 
tance  it  would  be  unreasonable  to  calculate  the  mean  square  error.  As 
a  measure  of  dispersion  we  choose  instead  the  differences  between  the 
highest  and  the  lowest  measured  values  (the  range).  The  greatest 
deviations  from  the  above  mean  values,  given  in  columns  2  and  4,  are 
thus  approximately  one-half  of  the  range  quoted  in  columns  3  and  5. 
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The  dispersion  in  r  varies  with  the  charge  weight 


TABLE  I 


MEASURED  TIME  FOR  SHOCK  FRONT  TO  REACH  r,  AND  TIME 
AT  MCE  THE  EXCESS  PRESSURE  VANISHES 

Values  corrected  to  1  kg  charge  weight,  0  deg  C  temperature,  and  0  per 
cent  humidii^r. 

Explosive  Charge:  TNT. 

r  =  0.054  m  kg 
o  * 


Q  =  0.25 

kg 

Q  =  1.00  kg 

B 

2 

3 

4 

5 

7 

8 

9 

10 

■ 

No. 

Dis- 

Dis- 

of 

r  -  r 

of 

r  "  ro 

persion 

t 

persion 

Obs . 

m*  0 

Obs. 

m 

TTTTn. 

ms 

//s 

Small  Distance: 

Wave  Front 

4 

.104 

1 

.021 

1 

4 

.200 

1 

.047 

2 

4 

.183 

1 

.043 

6 

4 

.300 

1 

.083 

2 

4 

.341 

2 

.100 

7 

4 

.400 

1 

.128 

6 

1 

.422 

- 

.133 

- 

5 

.446 

2 

.142 

30 

6 

.579 

1 

.213 

16 

5 

.546 

5 

.204 

19 

5 

.658 

4 

.264 

17 

5 

«  .746 

4 

.341 

21 

6 

.737 

4 

.330 

21 

4 

.846 

6 

.413 

15 

5 

.816 

4 

.389 

8 

4 

1.046 

2 

.511 

0 

6 

.897 

6 

.465 

27 

4 

1.246 

1 

•  C<i3 

20 

7 

1.056 

12 

.626 

31 

5 

1.215 

14 

.809 

48 

7 

1.374 

15 

1.018 

28 

5 

1.532 

10 

1.224 

21 

G 

1.690 

12 

1.491 

43 

5 

1.848 

13 

1.774 

24 

1 

2.001 

— 

2.062 

— 

Large  Distance: 

Wav 

3  Front 

2 

3.11 

10 

4.41 

90 

2 

3.95 

10 

6.54 

30 

2 

6.22 

20 

12.72 

0 

0 

5. 98 

30 

12.04 

10 

3 

9.44 

0 

21.79 

40 

2 

7.97 

10 

17.33 

40 

2 

12.57 

10 

30.68 

60 

2 

9.33 

20 

22.97 

50 

2 

15.74 

0 

39.81 

25 

2 

11.94 

0 

28.97 

40 

2 

18.35 

0 

49.22 

24 

Lar 

ge  Distance :  Time  at  which  Pos 

itive  Phase  vanishes 

for  Given  Distance  r 

2 

12.57 

10 

34.  90 

140 

2 

7.97 

10 

20.95 

450 

2 

15.74 

0 

44.37 

490 

2 

9.93 

20 

27.00 

90 

2 

18.85 

0 

53 . 65 

ISO 

2 

11.94 

0 

32.75 

70 

[sic;  the  dimensions  given  in  this  row  should  be  multiplied  by  kg 


9 


TABLE  II 


VARIOUS  AMBIENT  DENSITIES 
0  deg  C  temperature,  and  0  per 

Explosive  Charge:  TUT. 

Charge  weight:  0.25  kg. 

r  =  0.054  m  kg"1/3, 
o 


SHOCK  FRONT  TRAVEL  TIME  MEASURED  AT 

Values  corrected  to  1  kg  charge  weight, 
cent  humidity.  Density  »  in  kg  m"^ . 


r  -  r 

0 

Qo 

=  5.10 

JO 

o 

11 

3.81 

Qo  = 

2.55 

■O 

o 

II 

1.20 

M-  3 

of 

/* 

fxs 

% 

/^s 

at 

/« 

.15a 

50.6 

146 

45.3 

132 

40.8 

113 

34 . 6 

100 

.278 

— 

— 

99.9 

136 

— 

— 

73.5 

100 

.316 

131.7 

149 

121.3 

137 

105.9 

120 

88.2 

100 

.342 

147.6 

149 

136.4 

138 

— 

— 

98.8 

100 

.4  76 

249.0 

153 

229.1 

141 

202.0 

124 

162.6 

100 

.786 

534.0 

159 

534.6 

145 

— 

368.1 

100 

The  travel  tines  are  measured  by  one  shot  for  each  distance,  except  for 
p  =  2.55  where  two  shots  are  used,  ana  ^  =  1.29  where  the  tines  are 

calculated  by  the  formulae  for  small  distances. 

TABLE  III 

MEASURED  TRAVEL  TIMES  FOR  VARIOUS  EXPLOSIVE  CHARGES  COMPARED 
TO  THOSE  FOR  TNT  AT  ATMOSPHERIC  PRESSURE 

Mean  values  of  two  shots  at  each  distance. 
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1.05 

94 
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1.35 

94 
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95 
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1.65 
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100 

100 
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M.  v. 
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114 
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—  •  03 
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* 

** 


.  -1/5 
kg  in 


-  -1/z] 

sic;  m  .  kg  '  J 
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EVALUATION  OF  EXP3RIJSEHTAL  RESULTS 

The  measured  values  of  the  travel  times  and  t£  may  he  reprod\'.ced 
very  accurately  by  means  of  the  following  formulae; 

For  small  distances 

tl  =  Vr  -  r0)  +  Vr  -  r0)2  +  b2(r  "  ro'"; 

For  large  distances 

^2  =  c</(r  "  ro}  +  Cl/(r  "  ro)2 

where  cq  is  the  sonic  velocity,  i.e.,  33.3  m  sec.  ^  £sic;  331.3  m  sec 
The  b's  and  c's  are  determined  by  the  method  of  least  squares. 

Thus  the  following  values  are  obtained  for  a  spherical  TUT  charge; 

bQ  =  0.1614,  b1  =  0.3598,  bg  =  0.0390; 

c  =  331.3,  c.  =  375.7. 

o  1 

It  is  evident  from  Figures  4  and  5  that  the  measured  times  fall  on 
the  curve  computed  from  the  above  formula  quite  well,  and  that  values 
for  various  charge  v/eights  lie  on  the  same  curves  since  they  are  reduced 
in  conformance  with  the  laws  of  similitude  the  validity  of  which  is  thus 
confirmed  in  this  particular  «ase. 

•The  velocity  for  large  distances  is  calculated  fr  OTI1  Xi  -10  1710 jj  VuT  S  d. 
travel  times  by  the  method  of  Partlo  and  Service  (Reference  7). 

One  calculates  the  difference  A  between  the  travel  times  of  the 
sound  wave  (c  =  331.3  m  sec”l)  and  the  shock  wave  as  follows; 

A  =  (r  -  rQ)/c  -  t. 

Differentiating,  one  obtains 
d  A  /dr  =  l/c  -  l/.v 

whence 

w  =  c/  [l  -  c  (d  A  /dr  )J  . 
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■-r0"  •?  .«  .6  .8  1.0  1.3  1.1  nfcl-6 

Distance  from  Charge  Surface 

Travel  Times  of  the  Shock  Front  from  Different  Explosives 
Compared  to  TNT  at  Various  Distances  from  the  Charge 


FIG.  5 


Reciprocal  of  Distance  from  Charge  Surface 

Time  at  which  Positive  Phase  Ends  at  Various 
Distances  from  Spherical  TNT  Charges. 

+  £sio; 
*  laic; 
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The  derivatives  dA/dr  are  determined  by  smoothing 
:-r  -i  i  cn  1  Iv . 


iTi*3u.  jur  od  sTO-^VlSS 


3  “i.e  caicuxitev. 

•and  v/avs-lengths  for  a  spherical  TUT  charge 


D  i-  T>  l  -■•  '•*  .  a  ZTI w  o  y  iT  c *  a.  ^  j  a 

i.t  atmospheric  crscsur.D. 


The  measured  travel  tines  at  --arious  ambient  densities  or  “he 
air  are  evaluated  according  to  Tabic  lUby  comparison  with  the  corre¬ 
sponding  tines  for  T1IT  calculated  with  che  aid  of  the  formulae  mentioned 
before. 


TA5L3  IV 


CALCULATED  TRAVEL  TI15E3 
SPHERICAL  TUT 


,  FRCITT  VELOCITIES,  AID  m73  LEU 
CHARGES  AT  A'Z'OSPKERIC  PRESSURE 


GTRS  OF 
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Tlhen  r  -  rQ  ^  1.70,  t^  is  calculated  by  the  formula 

t±  =  0.1614  (r  -  r  )  +  0.3523  (r  -  r q)2  +  0.0390  (r 
(t^  in  milliseconds) 

r  -  rQ  ^  1.70,  tg  is  computed  by  the  formula 
1/t  =  331.3/(r  -  rQ)  +  373.7/(r  -  ro)2. 

(tg  in  sec 


'^ien  r  -  r q  ^  12,  one  has 

wL  =  331.3  +  160/r. 

(w^  in  m  sec 

The  wave  lengths  are  calculated  by  the  approximate  formula 


1 ,  =  c  t  . 
+  o  + 


These  dimensionless  times  are  shewn  in  Figure  6  as  functions  of 
(r  -  r  ).  For  a  given  density  q  the  values  apparently  lie  with 
sufficient  accuracy  on  a  straight  line.  Hence  one  may  extrapolate  with 
a-  satisfactory  degree  of  precision  to(r  -  rQ)  =  0.  •  Reciprocal  values  of 
the  travel  times  for  these  values  of  ^  gi-fe  the  ’ratio  of  the  shock 
front  velocity  w  at  the  density  under  consideration  to  the  velocity  at 
the  normal  density  of  1.29  kg  m”^  (see  Figure  7). 

It  is  apparent  from  Figure  8 
between  wQ  and  (>0  is  practically 
chooses  from  Table  IV  the  value  w 
formula 

wq  =  6190/(0.85  +  0.116  f}Q). 

(w  in  m  sec  M 

for  the  range  of  Q0  from  1.29  to  5.1.  The  present  experimental  data 
are  not  yet  sufficient  to  test  ?;hether  the  above  formula  holds  also  for 
the  negative  phase. 

The  travel  times  measured  at  atmospheric  pressure  for  various 
explosives  as  given  in  Table  III  are  also  evaluated  by  direct  comparison 
with  HIT  and  plotted  in  Figure  8,  which  shows  that  the  travel  times  and, 
therefore,  also  the  shock  front  velocities  of  each  explosive  relative  to 
TNT  are  not  constant  with  respect  to  the  distance  (r  -  r  ) . 


.t  the  relation 


r sic ;  Figure  7^  the 
linear  and  may  be  represented,  if  one 
o  =  6190  m  sec-^  for  %  o  =  1«29,  by  the 
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FIG.  6 


es 


Distance  from  Charge  Surface 


Travel  Time  at  Various  Densities  of  the  Surrounding  Atmosphere  at 
Various  Distances  from  SphericaL  TNT  Charges. 


FIG.  7 


?.  *  1  2  3  ♦  5  kg.a': 


Density  of  Ambient  Atmosphere 

Front  Velocity  at  the  Surface  of  the  Charge  in  Atmospheres  of  Various  Densities. 


FIG.  S 


Distance  from  Charge  Surface 

Travel  Time  of  the  Shock  Front  at  Various 
Distances  from  Spherical  TNT  Charges. 


SXAlHIIATIOi:  OF  RUDE 


ESRG’S  ATT  SEDOV’S  FORMULAS 


According  to  Rudenberg  the  shod:  front  velocity  v?  is  calculated 
from  the  energy  content  3  of  the  explosive  charge,  the  specific  density 
of  the  distributed  air,  and  the  ambient  density  by  the  formula 

vjq  =  (3  q'Q/  ^0)1//2. 

.  • 

Substituting  values  for  TUT,  this  formula  yields  a  value  of  v;Q 
which  is  about  one  third  of  the  measured  value.  For  a  given  explosive 
(S  (3’o  constant)  one  should  get 

w  CC  p 

o  no 

The  measured  values  for  TUT  do  not  obey  this  lam  at  all.  neither  does 
v/  change  proportionally  to  E^/~}  according  to  experiments.  Another 
method  to  approach  this  problem  has  been  presented  by  the  author  else¬ 
where  (Reference  3). 


Rudenberg  gives  for  the  relation  between  w  and  r 

-  °  ,-,1/2 
=[»  (  DA/'  1)  +  (/-  l)Zl/<8Jr$„r2)J  .. 

Therefore,  one  should  have  for  a  given  explosive  and  density 


where  a,  and  a„  are  constants.  This  means  that  w  “  should  be  a  linear 
12  o 

function  of  r“^. 


A  graphical  representation  of  these  values  according  to  Table  IV 
does  not  verify  this  law . 


aoeve, 


Sedov  gives  (Reference  2)  ^Appendix  aJ,  using  the  same  notation  as 

T1  “  <  •  «S  [»i«>  ®5«oA)lA]. 


lleither  of  the  above  dimensional  formulae  agrees  with  the  measurements. 

This  shows  in  addition  that  it  is  not  possible  to  express  I,  as  a  pro- 

j. 

duct  of  unique  functions,  in  the  form  T^  =  f^(  ()  q)  f  0  (h)  fj(R),  because 

f given  values  of  £>0and  3,  for  example,  the  ra*’  r>  n?  -'-r'.vel  -Mr-*** 

for  two  different  ambient  densities,  or  for  two  different  types  of 
explosives,  should  be  independent  of  rj  i.e.‘,  each  curve  in  both  Figures 
6  and  8  should  be  parallel  to  the  (r  -  r  )-axis,  whereas  actually  they 
are  not.  Therefore,  it  is  not  possible  to  find  a  general  distance  law 
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fj(R)  which  is  completely  independent  of  the  other  parameters  of  ex¬ 
plosion.  This  behaviour  is  explained  by  the  fact  that  changes  of  the 
quantities  describing  the  wave  front  are  affected  not  only  by  the  dis¬ 
tance  from  the  charge  but  also  by  the  pressure  distribution  behind  the 
front,  as  the  author  has  shown  elsewhere  (Reference  9).  This  pressure 
distribution  is  not  independent  of  the  properties  of  the  explosive 
material,  and,  therefore,  each  explosive  substance  must  have  its  own 
characteristic  pressure-distance  lav/. 

AC  IQICWL2D  GEMEtTT 

All  measuring  equipment,  with  the  exception  of  the  spiral  oscillo¬ 
graph  and  the  tubes  for  the  cathode-ray  oscillograph, was  designed  and 
built  by  the  Fysilcaliska  Forskingsavdelningen,  AB  Bofors.  The  measure¬ 
ments  were  made  with  great  care  by  Engineer  A.  Haggfcvist,  who, further¬ 
more,  proposed  and  developed  methods  to  provide  the  cathode-ray  oscillo¬ 
graph  with  a  FSTil  fuze  and  attended  to  the  numerous  important  details. 


SUMMARY 


Travel  times  of  shock  waves  emitted  from  bare  spherical  charges  of 
high  explosives  have  been  measured.  The  experimental  arrangements  are 
described.  Observed  values  may  be  reproduced  with  great  precision  by  . 
the  formulae  given  in  Figures  4  and  5. 

The  scaling  laws  have  been  verified  over  the  whole  range. 

The  relation  between  the  front  velocity  and  the  ambient  density  of 
the  surrounding  atmosphere  has  been  determined. 

The  front  velocities  of  shocks  from  PETIT  and  from  different  mixtures 
of  T13T  and  A1  have  been  determined  (Figure  8). 

The  resulting  data  do  not  verify  the  theories  of  Rudenberg  and 
Sedov. 


17 


REFERENCES 


1.  Rudeu’oerg,  R. 


2.  Sedov,  L.  I. 


3.  ITeibull,  W. 


4.  Wolff,  \7. 


5.  v.  Angerer,  S.,  and 
Ladenburg,  R. 


3.  Burlot,  E. 


7.  Partlo,  F.  I.,  and 
Service,  J.  H. 


liber  die  Fortpf lanzungs  go  sclnvindigkg  it  and 
Impulsstarke  von  Verdiclitungsstossen. 

Art.  Monath. ,  Vol.  113,  1916,  pp.  237-265, 
and  Vol.  114,  1913,  pp.  2S5-31S. 

[On  the  Velocity  of  Propagation  and  the 
Inpulse  of  Shock  Haves] 

Le  mouvenent  d’air  en  cas  d'une  forte  ex¬ 
plosion. 

C.  R.  Acad.  Sci.  DR33,  Vol.  52,  1946, 
pp.  17-20. 

[jhe  Motion  of  Air  in  a  Strong  Exp 1 os ion] 

Cm  den  mekaniska  verkan  av  brisar.ta.  s  prang - 
tinmen  vid  detonation  i  luft  och  i  vatten. 

Tidskr.  Kustart.,  ITo.  2,  1944 ,  pp.  3-35. 

[On  the  Mechanical  Work  of  Explosives  Detonated 
in  Air  and  in  Water] 

Uber  die  bei  Explosicnen  in  der  Luft  e  in  ge¬ 
ls  it  eten  VorgSnge. 

Ann.  Fhys.  Chen.,  Vol.  69,  1399,  p.  329. 

(On  the  Phenomena  Induced  by  Explosions  in  Ai^ 

Bxperimentelle  Beitrage  zur  Ausbreitung  des 
Schalles  in  der  Freien  AtmosphSre. 

Ann.  Fhys.,  Vol.  66,  1921,  p.  293. 

[Experimental  Contributions  to  the  Propagation 
of  Sound  in  Free  Ai£ 

Comnte  Rendu  Somnaire  des  Experiences  de  La 

m,m\m  ii  . .  ..  - .  ^  ..ii  - 

Courtine  et  des  Essais  qui  lea  ont  Precedees. 

Him.  Art.  Fran$.,  Vol.  4,  1925,  pp.  479-523. 

(Summary  Report  on  the  Experiments  at  La  Courtine 
and  the  Preceding  Experiments  at  other  Places] 

Instantaneous  speeds  in  air  of  explosion  re¬ 
ports  at  short  distances  from  the  source. 

Physics,  Vol.  S,  1925,  pp.  1-5. 


13 


8.  Weibull,  W. 


Boundary  Conditions  at  the  Surfa.ce  of 
Detonating  High  Explosives. 


Ark.  liat.  Astr.  Fys.,  Yol.  34  B,  No.  IS, 
1947. 

9t  Vfeibull,  W.  Ondes  planes  de  prossion  %.  front  discontinu. 

VTe  Gangr.  Intern,  Mi/c.  Appl.,  Paris,  1946. 
[Plane  Shock  Tfaves] 


19 


APPENDIX  A 


THE  M0TIQ1T  OF  AIR  IN  A  STRONG  EXPLOSION 
(Le  movement  d'air  en  cas  d’une  forts  explosion) 

*7 


L.  I.  Sedov 

Comptes  Rendus  (Doklady)  de  L’Academie  des  Sciences  de  l'HRSS 
1946,  Vol.  52,  No.  1,  pp. 17-20 

In  my  preceding  communications  (References  A-l,  A-2)  I  have  shown 
certain  new  important  families  of  solutions  of  non-linear  equations  in 
gas  dynamics  for  plane,  cylindrical  and  spherical  waves.  The  gas 
motions  corresponding  to  these  solutions  may  have  strong  discontinuties 
which  are  propagated  with  variable  velocity.  The  entropy,  therefore, 
assumes  different  values  for  different  gas  particles. 

It  is  of  interest  in  the  study  of  explosion  phenomena  to  determine 
the  law  of  variation  of  the  shock  wave  velocity  as  a  function  of  its 
position  with  respect  to  the  place  of  explosion  The  shock  wave  velocity 
and  the  initial  state  of  the  gas  determine  the  pressure  jumps,  the 
velocities  and  the  densities  of  the  gas  particles;  it  is  these  magnitudes 
that  should  define  the  mechanical  effect  of  the  explosive  action. 

Yfe  shall  now  consider  the  problem  of  the  propagation  of  an  ex¬ 
plosion  wave  under  the  following  assumptions: 

(1)  In  the  spherical  case  one  neglects  the  dimensions  of  the 
charge  and  we  consider  it  as  a  point;  the  energy  liberated  by  the  ex¬ 
plosion  is  regarded  as  finite. 

(2)  At  the  initial  moment  (t  =  0)  air  is  at  rest.  At  the 
centra  of  explosion  energy  is  released  instantaneously  and  the  air 
traversed  by  the  shock  front  is  disturbed. 

(3)  Atmospheric  pressure  is  neglected  relative  to  the  shock 
front  pressure. 

The  last  assumption  shows  that  the  results  of  the  solution  do  not 
depend  on  the  initial  pressure  of  the  gas,  but  may  depend  on  the  initial 
density  0  which  characterizes  the  inertia  of  the  gas  particles. 

The  assumed  ideal  regime  does  not  take  into  account  factors  7/hich 
influence  the  phenomena  taking  place  in  the  immediate  vicinity  of  the 
charge  (the  finite  dimension  of  the  charge,  the  mass  of  the  charge,  etc.). 
At  a  large  distance  from  the  point  of  explosion  the  shock  wave  becomes 
a  sound  wave;  in  this  process  the  atmospheric  pressure  is  important  and 
now  the  third  assumption  is  inadmissible  for  the  description  of  the 
phenomena  at  large  distances  from  the  charge. 
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The  above  assumptions  may  be  considered  reasonable  where  the 
shock  is  still  strong  but  where  the  distance  is  large  relative  to  the 
charge  radius.  Under  such  assumptions  one  may  choose  the  following 
independent  determinative  dimensional  parameters  for  the  problem  of 
adiabatic  motion  of  an  ideal  gas  : 

the  time  t; 

the  coordinate  r  representing  the  distance  from  the  point  of) 
explosion; 

the  atmospheric  density 

a  constant  S,  proportional  to  the  energy  liberated  by  the 
explosion. 

One  may  modify  the  equation  of  motion  and  the  equation  of  state  of 
the  gas  to  the  system  of  independent  quantities  determining  the  dimension. 

Now  one  may  introduce  different  important  physical  constants  the 
dimensions  of  which  are  expressible  by  the  dimensions  of  S  and  of 

Having  defined  the  system  of  characteristic  parameters  in  this 
manner  by  means  of  dimension  theory,  it  is  simple  to  obtain  the  formulae 
which  yield  the  laws  of  variation  of  the  shock  front  velocity  c  and  of 
the  shock  front  distance  r*. 

In  the  spherical  case  one  has 

r*  =  (E/£q)1//5  t2/5,  c  =  dr*/dt  =  (2/5)  (E/^)1/2  .(r*)"3/2;  (l) 
similarly,  in  the  cylindrical  case 

r*  =  (E/<oo)1//4  t1/2,  c  =  (1/2)(E/^q)1/2  (r*)"1;  (2) 

and  in  the  plane  case 

r*  =  (E/ <Oo)1//3  t2/3,  c  =  (2/3)(E/(Oo)1/2  (r*)"1/2.  (3) 

It  follows  from  the  formulae  (l),  (2),  and  (3)  that  for  differently 
shaped  charges  one  obtains  different  laws  of  variation  of  the  shock  wave 
velocity  with  distance  from  the  point  of  explosion. 

At  points  behind  the  shock  front  we  have 

v=(r/t)V(A),  ^=^0R(A),  P  =  (O  Q(r/t)2  P(A  ).  (4) 

In  the  spherical  case  the  abstract  parameter  A  is  defined  by  the 
formula 


A  =  (E/£Q)(t2/r5). 

By  (4)  the  problem  of  determining  the  motion  of  a  gas  nay  be  reduced 
to  one  of  ordinary  differential  equations.  If  the  gas  is  ideal  and  the 
process  adiabatic  in  the  domain  of  continuous  motion,  the  3.bove  equations 
become  in  the  case  of  spherical  symmetry  (Reference  A-2) 
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d2  _  [2(7-1)  +'3(f-l)j  (7-  |)  -  (f-l)Y(7-l)(Y-|)  -  [2(V-1)  + 

dY  "  |>  (7-  |)  +  (57  “  3  Y)  ZJ  (6) 

d  in  A  1/5  _  (Y  "  ?)  “  Z _  f7, 

d  V  7(7-1)  (7  -  |)  +  (-gy  -  3  7)  z 


d  111  -  |)  R  .  JY 

d  In  A1//5  7  - 


(3) 


where  z  =  y  ?/l  and 
heat  capacities. 

It  nay  be  shewn  that  the  system  of  equations  (6),  (7)  and  (8) 
possesses  an  integral 

Y2  (7  -  |)  a  +  yrr  ?  (V  -  |)  +  2  FY  =  aA  (9) 

where  A  is  a  certain  constant. 


CQ/C 


the  ratio  of  the  isobaric  and  isochoric 


To  the*  state  of  unperturbed  air  corresponds  a  particular  solution 
of  the  equations  (6),  (7)  and  (8)  with  z  =  0;  7  =  0;  £3  =  ^  ;  R  =  1  and 
A  arbitrary.  Across  the  shock  front  the  state  changes  abruptly  and  one 
has  (Reference  A-2) 

7  =  0.8/(  f  +  1);  z  =  0.32  £  (  y  -  l)/(  Y  + 

R  =  (  ^  +  l)/(  %  -  l)j  A  =  const. 

The  constant  A  vanishes  for  the  state 
(10).  In  this  case  the  integral  (9)  leads 
equation  (6) 

2  _  (  y  -  1) (0.4  -  7)  y2  _ 

2  (v  -  §  ? ) 

To  the  centre  of  symmetry  corresponds 
equation  (6),  namely 

7  =  2/7  and  z  =  00. 

Using  (ll),  we  may  express  A  and  R  in 
7  by  starting  from  the  equations  (7)  and  (8).  In  this  manner  we  obtain  a 
simple  and  complete  solution  of  the  problem. 

Figures  A-l,  A-2  and  A-3  show  the  laws  of  variation  of  the  velocity, 
the  density  and  the  pressure  of  air  behind  a  shock  wave  {  Y  ~  l»-)j  ‘fcbe 
constant  S  is  related  to  E, ,  the  energy  liberated  by  the  exnlosicn,  by 
E  =  1.175  E1. 


corresponding  to  the  formulae 
to  the  following  solution  of 


(11) 

a  particular  solution  of 

finite  form  as  functions  of 


22 


At  the  centre  of  the  explosion  the  velocity  and  the  density  are 
zero;  the  pressure  approaches  zero  as  the  tine  increases.  This  shows 
that  when  a  wave  becomes  weak  an  inverse  motion  of  the  gas  toward  the 
centre  of  the  explosion  must  take  place. 

By  means  of  the  equations  (6),  (7)  and  (8)  one  may  obt3.in  the 
solution  of  (4)  ivith  several  or  even  infinitely  many  shocks. 
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fr£7f* 


APPENDIX  B 


COMPARISON  WITH  C0LTIL3D  DATA 

Data  on  spherical  Pentolite  compiled  at  Aberdeen  have  been  corrected 
from  an  assumed  temperature  of  20  deg  C  to  0  deg  C  and  are  plotted  on 
Figure  B  as  a  smooth  curve. 

* 

The  corresponding  part  of  Weibull's  data  on  spherical  TNT  for  Q  =  0.25 
kg  as  given  in  Table  I  is  plotted  on  the  same  figure. 
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